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FOREWORD 


At  the  present  time,  thousands  and  thousands  of  people  around  the  world  deal  with  ice, 
snow  and  permaftost  They  are  scientists,  educators,  engineers,  navigators,  meteorologists 
and  others.  While  a  anall  fraction  of  these  people  cont^ute  to  the  knowledge  base  in  ice 
physics,  all  of  them  use  knowledge  from  it  frequently.  Moreover,  successful  applied  research 
is  based  upon  fundamental  science — one  more  reason  for  ice  specialists  to  have  a  textbook 
on  ice  physics  on  their  desks. 

The  first  modem  ice  physics  text  was  Fletcher's  book  on  The  Chemical  Physics  cflce  (1970). 
Fletcher's  book  is  in  typical  textbook  fomnat:  it  is  teasoruibly  brief  and  easy  to  understand. 
He  toudied  rm  a  few  of  the  most  important  topics,  but  not  all  of  them. 

The  most  crunprehensive  book  on  ice  physics  to  date  was  published  by  Hobbs  in  1974. 
Hobbs  considered  almost  all  of  fiie  basic  aspects  of  ice  as  understood  at  drat  time.  Moreover, 
he  described  and  compared  several  (sometimes  opposing)  viewpoints.  This  fundamental 
and  rather  large  (837 pages)  book  is  conunoitly  known  as  the  'Tee  Bible"  by  specialists  in  the 
field.  In  1974  and  1975,  two  CRREL  monographs  on  ice  were  produced  by  John  Glen.  These 
were  briefly  and  clearly  written  and  reviewed  almost  all  ice-related  subjects.  This  work  was 
(and  in  some  respects  still  is)  a  magnificent  introduction  to  ke. 

Finally,  in  1981  Maeno  wrote  a  simple,  popular  book  for  die  e]q>ress  purpose  of  attracting 
people's  attention  to  the  subject 

During  the  past  20  years,  a  significant  amount  of  new  experimental  and  theoretical  work 
has  appeared,  dramatically  changing  our  views  on  ke  physics.  As  a  result  we  are  now  able 
to  formulate  physical  laws  using  more  simple  and  direct  methods.  We  have  found  some  of 
the  physical  models  used  in  dte  past  to  be  ctonpletely  wnmg.  The  physics  of  ice  is  a  much 
better  developed  subject  than  it  was  20  years  ago. 

For  die  above  reasrms,  we  feel  the  time  is  ripe  for  a  contemporary  book  on  ice  physics, 
incorporating  the  known  and  proven  with  almost  20  years'  worth  of  material  not  coveted  by 
previous  works. 

We  have  tried  to  prepare  a  "readable"  book,  and  not  one  that  requires  the  reader  to  be  a 
vmiquely  educated  person.  It  is  our  intent  to  present  die  material  in  such  a  way  that  any 
reader  attracted  by  ^  tide  Ice  Physics  will  be  abte  to  comptdiend  it  This  is  quite  difficult  for 
a  book  dedicated,  not  to  a  particular  field  of  knowledge,  but  to  a  specific  material.  Indeed, 
for  ice  it  means  we  have  to  consider  a  wide  variety  of  subjects,  including  quantum  chemistry, 
solid  state  physics,  the  theory  of  elasticity,  ionk  conductivity,  synchrotron  x-ray  topogra¬ 
phy,  crystal  growth,  the  physics  of  surfaces  and  more. 

The  primary  goal  is  to  produce  as  simple  a  book  as  possible  widiout  sacrificing  scientific 
accuracy.  Experimental  facts,  physical  ideas  and  dieories  will  be  strongly  organized  and 
bound  togedier  cohesively.  The  reader  will  be  introduced  to  a  wide  variety  of  material  on  a 
step-by-step  basis.  Then  die  picture  will  be  wh(de. 

To  accelmte  publication,  diis  book  will  appear  first  in  the  form  of  a  series  of  joint  CRREL- 
Dartmouth  tepotts,  later  to  be  published  in  CRREL's  Monograph  series,  on: 

1.  The  structure  of  ordinary  ice 
Part  I:  Ideal  structure  of  ice.  Ice  crystal  lattice 
Part  Ik  Defects  in  ke 
Volume  1:  Point  defects 
Volume  2:  Didocatiens  and  plane  defects 


2.  Electrical  properties  of  ice 

Part  I:  Conductivity  and  dielectric  permittivity  of  ice 
Part  II:  Advanced  topics  and  new  physical  phenomena 

3.  Optical  properties 

4.  Electro-optical  effects  in  ice 

5.  Thermal  properties 

6.  Mechanical  properties  of  ice.  Elasticity  and  anelastic  relaxation.  Plastic  properties. 
Frachue  of  ke 

7.  Elect^omechanic^ll  effects  in  ice 

8.  Surface  of  ice 

9.  Other  forms  of  ice  and  their  properties 

10.  Ice  in  space 

11.  Ice  research  laboratories 

The  reports  will  be  prepared  in  a  sequence  convenient  to  tire  autirors.  The  present  is  the 
fourth  in  the  series. 
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NOMENCLATURE 

This  Nomettclature  section  incorporates  that  of 
Volume  1,  Point  D^ts. 

a  lattice  constant 
b  Burgers  vector 
b  magnitude  of  Burgers  vector 
c  lattice  constant 
D  diffusion  coefficient 
Du  diffusion  coefficient  of  hydrogen  in  ice 
Dj  diffusion  coefficient  of  interstitials 
O,  self-diffusion  coefficient 
Dw  (wi^)exp(^j 
Oy  diffusion  coeffidoit  of  vacancies 
D  electric  displacement  vector 
e  proton  charge 
c  relative  dielectric  permittivity 
e  strain 

E^b  energy  of  creation  of  a  pair  of  Bjerrum 
defe^ 

Eai  energy  of  creation  of  an  ion  pair 
Em  activation  energy  of  static  conductivity  o, 
Em.  activaticm  energy  of  high-frequency  con¬ 
ductivity  o.. 

Ef  energy  of  a  defect  creation 
Oi  defects'  electric  charge  (t  =  1, 2, 3, 4) 

Ejf  formation  energy  of  interstitials 
E](  Idnk  formation  oiergy 
En^  activation  energies  of  protonic  defects' 
motkm(i  =  l,2,3,4) 

Cg  dielectric  permittivity  of  vacuum 
£,  activation  energy  of  self-diffusion 
e,  static  dielectric  permittivity  (ox  <  cod) 

Eoi  activaticHi  energy  of  partial  conductivity 
of  i-type  defects  (i  =1, 2, 3, 4) 

Ey  activationenergyofDebyerelaxationtime 
E  activaticHi  energy 

E..  hig^ftequetKy  (ca  <  <  ooq)  dielectric  per¬ 
mittivity 

e  strainrate 
F  freeenergy 
9  s  S^SkaTfoo 
/,v  frequencies 
I  fault  vector 


free  energy  of  defect  formation 
Fit  free  energy  of  kink  formation 
Fm  free  energy  of  kink  motion 

force  acting  on  t-type  of  defects 
g  diffraction  vector 

activation  volume  of  formation  of  pro- 
toruc  defects 

Ymi  activation  volume  of  protonic  defects' 
motion 

A  separation  of  the  Peierls  troughs 
Hi  =  1,-1,— 1,1  for  i  =  1, 2, 3, 4 
Ji  flux  density  of  defects  (i  =  1, 2, 3, 4) 
displacement  current 
drift  current 

J,J  electric  current  density 
Boltznumn  omstant 

1  unit  vector  parallel  to  a  dislocation  line 
p  mobility 

m  mobility  of  i-type  of  defects 
n  crmcentration  (in  m'®) 

Nj  concentration  of  dislocations 
Hq  D-defectamcentration 
nHzO  concentration  of  water  molecules  in  ice 
Hi  concentration  ofdefects(i  =  1,2, 3, 4) 

»!(  number  of  kinks  per  unit  lengffi 
L-defect  concentration 
rty  concentrations  of  vacancies 
P  pressure 

electric  charge  of  carriers 
Q  rate  of  heat  generation 
r  distance 

Too  oxygen-oxygen  distance  in  ice  lattice 
(2.76A) 

S  entropy 
<7  conductivity 
a  normal  stress 

Ob  Bjerrum  defect  conductivity — 0B=(%-f(T4 
Sf.  configurational  entropy 
Sf  vibrational  entropy 

partial  conductivity  of  t  charge  carrier 
Oion  ionic  conductivity — =  <^  +  (% 
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S|(  kink  entropy 

a,  static  or  low-frequency  conductivity  (oa  < 
<<0d) 

high-frequency  conductivity  (eo  >  >  (Od) 

T  temperature 
X  ^ear  stress 
t  time 

tt,  mean  time  to  reorient  a  hydrogen  bond 
Xq  =  Debye  relaxation  time 
lifetime  of  charge  carriers 
U  internal  energy 

Ujg^  activationenergyofmotionofinterstitials 
activation  energy  of  iomc  motion 


dislocation  velocity 
»d  drift  velocity 

activation  volume  of  formation  of  inter¬ 
stitials 

Vim  activation  volume  of  motion  of  interstitials 
kink  velocity 

Vmoi  moleciilar  volume  in  ice  (3.3  x  10“^  cm^) 
W  number  of  configurations  of  a  system 
to  circular  frequency 
(1^  Debye  frequency 
n  configurational  vector 
*  * 

y  =  mq/ nH2o 


Structure  of  Ordinary  Ice  Ih 
Part  II:  Defects  in  Ice 
Volume  2:  Dislocations  and  Plane  Defects 

VICTOR  F.  PETRENKO  AND  ROBERT  W.  WHITWORTH 


INTRODUCTION 

Dislocations  ate  line  defects  in  crystals  and  are 
most  easily  visualized  in  terms  of  fte  slip  of  one 
plane  of  atoms  over  another  by  a  single  lattice 
spacing,  as  illustrated  in  Figure  1.  In  the  interme- 
(^te  stage  shown  in  this  figure,  the  line  BC  is  the 
boundary  betweoi  the  slipped  region  ABC  and 
die  remainder  of  the  slipplane.  This  boimdary  line 
is  ^dislocation  line,  ot dislocation  for  short.  The  pres¬ 
ence  of  the  dislocation  results  in  an  elastic  distor¬ 
tion  diroughout  the  whole  crystal,  but  the  disrup¬ 
tion  of  the  lattice  is  greatest  along  the  core  BC.  The 
lattice  vector  by  whidr  slip  has  occurred  on  the 
plane  ABC  is  known  as  the  Burgers  vector  b  of  the 


dislocation.  Where  the  dislocation  lies  parallel  tob 
as  at  B,  it  is  said  to  have  screw  orientation,  and 
where  it  is  perpendicular  to  b  as  at  C,  it  has  edge 
orientation;  in  general,  a  dislocation  will  have  both 
screw  and  edge  components.  An  edge  dislocation 
has  dangling  bonds  at  its  core,  as  seen  in  the 
enlarged  region  around  C  shown  in  Figure  1,  and 
can  be  thought  of  alternatively  as  the  line  at  which 
an  added  half-plane  of  atoms  terminates  within 
the  crystal.  There  are  no  dangling  bonds  on  a 
screw  dislocation,  but  the  lattice  is  distorted  in 
such  a  way  that  a  path  passing  once  round  the  core 
advances  by  one  Burgers  vector. 

From  the  way  in  which  we  have  introduced 
them,  it  is  clear  ^at  the  motion  of  dislocations  on 


Figure  1.  Production  of  slip  in  a  crystalline 
material  by  theglideofadi^ocation  BC  across  the 
slip  plane.  The  amount  of  slip  is  given  by  the 
Burgers  vector  b.  The  magnified  portions  show 
thearrangementofatomsaround  thescrewcom- 
ponent  at  Band  the  edge  component  ate. 


the  slip  plane  is  associated  with  plastic  deformation 
of  the  crystal.  Such  motion  of  the  dislocation  is 
known  as  glide,  but  it  represents  only  one  aspect  of 
die  properties  of  dislocations.  Elislocations  may  be 
incorporated  into  a  crystal  as  it  is  grown,  aff^ting 
the  topology  of  the  whole  lattice,  but  with  little 
eff^  on  plastic  deformation.  The  core  of  a  disloca¬ 
tion  may  be  displaced  perpendicular  to  the  glide 
plane  by  the  adding  of  atoms  to  or  removing  of 
atoms  f^m  tiie  end  of  the  extra  half-plane;  this  is 
known  as  climb. 

The  properties  of  dislocations  are  described  in 
many  l^le  and  articles.  For  a  very  clear  discus¬ 
sion  of  the  basic  geometry,  we  recommend  Read 
(1953),  and  for  a  comprehensive  discussion  includ¬ 
ing  many  aspects  relevant  to  ice,  Hirth  and  Lothe 
(1982).  A  beautifully  illustrated  book  describing 
dislocations  in  ice  is  that  edited  by  Higashi  (1988). 

DISLOCATIONS  IN  THE  ICE  STRUCTURE 
Basal  dislocations 

It  is  now  well  established  that  crystals  of  ice 
deform  by  slip  on  die  basal  plane  (0001)  (Glen  and 
Perutz  1954),  and  that  macroscopic  slip  on  any 
other  plane  is  difficult  (Higashi  1%9,  Duval  et  al. 
1983).  The  Burgers  vectors  for  slip  on  the  basal 
pl^  are  the  three  lattice  vectors  of  the  form  (a/3) 
<2110>,  but  in  macroscopic  experiments  slip  can 
occur  in  any  direction  by  a  combination  of  di^oca- 
tions  with  tiiese  tiuee  vectors  (Kamb  1%1).  We  will, 
therefore,  consider  initially  only  basal  dislocations 
of  this  type.  Because  of  ttie  hexagonal  symmetry, 
the  simplest  dislocatimis  ate  those  that  lie  parallel 
or  at  60°  to  the  Burgers  vector;  tiiese  ate  screw  and 
60°  dislocations  respectively. 


Figure!.  Projectionofthestructurecfice(mthe{121Q} 
plane,  showing  the  himl planes  cfthe shuffleset  and  of 
theglideset;landllared^perentkindsof{1010)planes. 
Only  oxygen  atoms  are  shown. 


J^igure  2  shows  the  ice  structure  projected  on  a 
(1210)  plane  with  the  bcisal  plane  horizontal.  Slip 
can,  in  principle,  take  place  between  two  kinds  of 
planes  called  the  glide  set  and  the  shuffle  set.  The 
planes  of  the  shuffle  set  are  more  widely  spaced 
and  were  originally  thought  to  be  the  natur^  slip 
planes,  but  the  planes  of  the  glide  set  tit  over  one 
another  in  a  way  that  resembles  the  packing  of 
close-packed  metals.  O;  such  planes,  a  dislocation 
may  lower  its  energy  by  dissociating  into  two 
partial  dislocations,  separated  by  a  stacking  fault 
(see  tile  Stacking  Faults  section).  This  dissociation 
is  illustrated  in  Figure  3a,  in  which  the  open  and 
shaded  circles  represent  oxygen  atoms  on  adjacent 
basal  planes  of  the  glide  set;  the  screw  dislocation 
in  Figure  3b  is  shown  dissociating  into  two  partial 
dislocations  j>f  Burgers  vectors  bi  and  \>2  of  tiie 
type  (a/3)  <1100>.  Slip  by  bi  corresponds  to  mo¬ 
tion  of  the  shaded  atoms  in  Hgure  3a  from  B  to  C, 
creating  the  stacking  fault,  and  slip  by  the  furtiier 
amount  b2  removes  the  fault,  resulting  in  a  net  slip 
by  the  full  Burgers  vector  b  =  bi  +  b2.  Both  tiie 
partial  dislocations  associated  with  a  screw  dislo¬ 
cation  have  a  30°  character.  Figure  3c  shows  the 
dissociation  of  a  60°  dislocation  into  an  edge  and  a 
30°  partial  dislocation.  The  dissociation  leads  to  a 
reduction  in  the  elastic  strain  energy  of  the  dislo¬ 
cation,  and  this  reduction  is  balanced  by  the  ener¬ 
gy  required  to  create  the  stacking  fault.  From  esti¬ 
mates  of  the  stacking  fault  energy,  Fukuda  et  al. 
(1987)  have  estimated  that  in  ice  screw  disloca¬ 
tions  cm  planes  of  the  glide  set  would  dissociate 
into  partial  dislocations  separated  by  about20nm. 

In  many  semiccmductors  with  structures  relat¬ 
ed  to  that  of  ice,  the  disscxiation  of  dislcxations  on 
tiie  basal  plane  has  been  cibserved  directly  in  the 
electron  microscope,  andslip  istiierefore  assumed 
to  (xcur  on  planes  of  the  glide  set  (see  George  and 
Rabier  1987).  There  is  no  sudi  direct  evidence  for 
ice;  and,  as  the  bonding  is  quite  different,  we  can¬ 
not  presume  that  the  same  ccmclusion  applies. 
However,  as  we  will  see,  indirect  evidence  sug¬ 
gests  that  disl(x;ations  do  in  fact  glide  on  planes  of 
foe  glide  set. 

If  a  straight  dislocation  lying  on  its  glide  plane 
moves  forward  by  one  Burgers  vector  over  part  of 
its  lengtii,  the  step  so  produced  is  caUed  a  Jb'nk.  An 
example  for  an  edge  dislcxation  is  shown  inF^ure 
4.  The  most  elementary  step  in  the  process  of  glide 
is  for  a  kink  to  move  cdcmg  the  dislocaticm  by  cme 
lattice  spacing,  as  shownby  thebroken  line.  Rgure 
4  also  shows  a;og  at  which  the  dislcxation  makes 
a  step  from  one  glide  plane  to  another.  For  a  dislo¬ 
cation  with  an  edge  component,  the  formation  of  a 
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Figure  3,  Dissociation  into  tioo  partial  dislocations.  Screw  and  60°  dislocations  on  planes  of  the  glide  set  can 
dissociate  into  partial  dislocations  with  the  Burgers  vectors  bj  and  bz  separated  ty  a  staddng  fault,  shown 
shaded. 

jog  or  its  motion  along  the  dislocation  requires  that 
atoms  oe  added  to  or  removed  ht>m  the  end  of  the 
extra  half-plane,  and  this  is  what  is  involved  in 
climb.  In  this  process,  the  dislocation  acts  as  a 
source  or  sink  of  vacancies  or  interstitials.  In  the  ice 
structure,  the  simplest  jog  would  take  the  disloca¬ 
tion  from  a  glide  to  a  shuffle  plane;  a  full  jog  must 
span  two  layers  of  atoms  (see  Fig.  2). 

Kinks  and  jogs  are  in  general  quite  different  in 
character.  However,  a  pure  screw  dislocation  has 
the  property  that  it  can,  in  principle,  glide  on  any 
plane  containing  its  Burgers  vector,  and  in  this 
case  whether  a  particular  step  behaves  as  a  kink  or 
a  jog  depends  upon  which  plane  is  being  consid¬ 
ered  to  be  the  glide  plane. 

Figures.  Edgedislocationcontainingakinkon  its  glide  Figures  5a  andb  ^ow  the  arrangement  of  mol- 

plane  and  a  jog  at  ivhich  it  nukes  a  step  from  one  glide  ecules  on  the  two  sides  of  the  basal  glide  plane  at 

plane  to  another.  The  broken  line  shows  an  elementary  an  edge  and  a  30°  partial  dislocation  in  the  ice 

step  in  the  motion  of  the  kink,  structure  (Whitworth  1980).  Open  circles  repre- 
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- ►[1120]  - ►[1120] 

a.  Dishcatum  con  for  edge  dislocations  vHth  dangling  b.  Dislocation  con  toiOt  30“  partial  dislocations  and 

bonds.  dangling  bonds. 


c.  Possible  way  in  which  the  con  in  Figun  5a  can  be  d.  Possible  way  in  which  the  con  in  Figun  5b  can  be 
reconstructed  to  link  up  the  dangling  Imds.  reconstructed  to  link  up  the  dangling  bonds. 

Figure  5.  Positions  cf  molecules  idHfoe  and  below  the  (0001)  glide  plane  for  partial  dislocations  in  ice. 


sent  molecules  above  the  glide  plane  and  shaded  Nonbasal  dislocations 

didesiepresentdiosebelow.Theextrahalf-plane  Disloc^ons  with  the  basal  Burgers  vector 

is  above  the  g^ide  plane,  and  the  dislocations  in  (a/3)  <2110>  can  in  principle  glide  on  the  non- 

both  of  die  diagrams  include  kinks.  There  are  bas^  planes  Jhat  include  diis  vector,  sudi  as 

dangling  bonds  on  the  molecules  above  the  glide  (0110)  or  {0111),andinthe£)isloaitionsAssodated 

plane  at  the  edge  dislocation,  and  <Hi  both  sides  of  Plastic  Damnation  section,  we  will  report 

the  glide  plane  for  die  30°  dislocatirm.  With  some  observations  of  edge  dislocations  diat  glide  on 

changes  to  bond  lengdis  and  bond  angles,  and  suchplanes-Thereisagainmorethanonepossible 
with  some  local  elastic  distortion,  it  is  possible  to  typeofglideplane,suchasd)esetsof  (lOlO)planes 
link  up  diese  dangling  bonds,  except  at  the  kinks,  marked  I  and  n  in  Rgure  2,  but  dteie  is  no  evidence 

as  shown  in  Figures  5c  and  d.  This  is  known  as  to  show  whkh  of  these  is  operative, 

reconstructibn,  whkhisgenerallybelievedtooccur  biadditimiwemayexpectdislocationsinwhich 

in  semiomductors.  Calculations  by  Heggie  et  al.  the  Burgers  vector  bias  a  [0001]  component.  Such 
(1992)  suggest  dmt  it  probably  happens  in  ice,  too,  dislocati<»is  have  been  observed  in  ice,  usually  in 

aldiou^  there  is  more  than  one  possibility  about  the  form  of  loops  lying  in  the  (0001)  plane.  Such 

how  the  recxHistruction  might  occur  for  die  30°  loops  are  called  prismatic  and  are  formed  by  die 
dislocaticHi.  condensation  of  point  defects.  There  is  no  evi- 


dence  that  these  dislocations  can  glide  or  that  ice 
can  be  deformed  plastically  by  slip  in  the  [0001] 
direction. 

DIRECT  OBSERVATION  OF 
DISLOCATIONS 

General 

Since  the  earliestexperiments  of  Webb  and  Hayes 
(1967),  X-ray  topography  hasbeen  extensively 
for  the  study  of  dislocations  in  ice.  It  is  undoubtedly 
dw  most  suitable  teduiique  for  use  vvidi  this  mate¬ 
rial,  and,  of  all  materials,  ice  has  been  the  one  most 
fruitfully  studied  by  this  method.  We  will  first  try 
to  explain  why  diis  is  so. 

For  other  materials  the  most  powerful  technique 
for  observing  dislocations  is  usually  the  transmis¬ 
sion  electnm  microscope.  However,  in  the  case  of 
ice,  the  preparaticxi  and  handling  of  suitable  thin 
specimens  present  enormous  difficulties,  and,  even 
when  prepared,  the  specimens  have  a  very  limited 
life  in  ttie  electron  beam.  Although  there  lui  ve  been 
reports  of  such  experiments  (Unwin  aiul  Muguru- 
ma  1972,  Falls  et  al.  1983),  no  significant  informa¬ 
tion  has  been  obtained  in  this  way.  It  is  possible  to 
produce  etch  pits  on  ice  (Higuchi  1958)  that  can  be 
useful  in  orienting  cr^tals,  but  these  do  not  nor¬ 
mally  correspond  to  the  points  of  emergence  of 
dislocations.  Particular  kinds  of  etch  pits  and  other 
etch  features  do  appear  to  be  related  to  dislocations 
(e.g.,  Mugururtui  and  Higashi  1963,  Sinha  1977), 
but  toe  information  obtained  from  such  experi¬ 
ments  is  very  limited  compared  with  that  from 
ottier  crystals,  such  as  UF. 

In  contrast.  X-ray  topography  has  been  used  to 
reveal  dislocations  in  toe  interior  of  crystals  that  are 
a  few  millimeters  thick  and  to  observe  their  motion 
during  annealing  and  while  under  stress.  This  is 
possible  because  ioe,havingalowmolecular  weight, 
is  sufficiently  transparent  to  X- 
rays  of  wavelengto  less  than 
al^t0.9A,  and  ice  crystals  can 
be  grown  with  a  low  enough 
dislocation  denaty,tiiatmdivid- 
ual  dislocations  can  be  distin¬ 
guished  in  suchlargeq>ecimens. 

A  furttier  reason  why  ice  is  suit¬ 
able  for  dyruurtk  experiments 
on  dislocatkms  is  that  titey  can 
be  moved  slowly  tmder  stress, 
whereas  in  many  materials  dis¬ 
locations  move  suddenly  by 
large  distances  once  some  criti¬ 
cal  stress  is  reached. 


X-ray  topography  technique 

The  form  of  X-ray  topography  that  is  the  sim¬ 
plest  and  the  most  easily  explained  has  only  be¬ 
come  possible  with  the  availability  of  intense, 
highly  collimated  beams  of  "white"  X-radiation 
from  a  synchrotron  source.  If  such  a  beam  falls  on 
a  single  crystal,  as  shown  in  Figrue  6,  it  produces 
Laue  diffraction  sp>ots,  in  which  there  is  a  one-to- 
one  correspondence  between  a  position  on  toe 
spot  and  toe  position  at  which  toe  X-rays  passed 
through  toe  aystal.  Lcxal  misorientations  of  the 
lattice  within  toe  crystal,  such  as  cxxur  at  dislcxa- 
tions,  change  the  difiraction  conditions  slightly 
and  result  in  contrast  in  the  Laue  spot.  E2ich  spot  is 
toerefore  an  image  of  the  crystal  in  which  disloca¬ 
tions  are  visible,  and  these  images  are  called  topo¬ 
graphs.  The  topographs  can  be  recorded  on  high- 
resolution  photographic  film  or  plates,  and  ob¬ 
served  in  real  time  at  lower  resolution  with  an  X- 
ray  sensitive  TV  camera. 

The  imaging  conditions  depend  on  the  diffrac¬ 
tion  vector  g,  so  that  a  dislcxation  of  Burgers  vector 
b  willnot  appear  in  a  topograph  for  whi^both  g-b 
= 0  and  g-(b  X  I) = 0,  wh^  1  is  a  unit  vector  parallel 
to  the  dislcxration  line.  This  nrakes  it  possible  to 
identify  ttw  character  of  each  dislocation. 

The  resolution  depends  on  toe  degree  of  colli- 
matim  of  the  incident  beam,  which  can  only  be 
satisfactorily  achieved  with  a  synchrotron.  With¬ 
out  this  collimation,  we  have  to  use  monochro¬ 
matic  radiaticm,  and,  with  a  ccmventional  X-ray 
source,  toe  beam  divergence  for  a  reasonable  in¬ 
tensity  is  then  too  high  for  Bragg's  law  to  be  satis¬ 
fied  over  the  whole  crystal  at  once.  Provision  has 
then  to  be  made  for  scanning  toe  crystal  across  toe 
beam,  and  toe  most  commc^y  used  arrangement 
is  tite  Lang  (1959)  camera,  which  has  been  used 
extensively  in  all  but  toe  most  recent  work  on  ice 
(see  review  by  Higashi  1988).  Because  of  the  scan- 


ning,  d)e  exposure  time,  which  is  typically  a  few 
seconds  using  a  syndm>tron,  becomes  minutes  or 
hours,  depofiding  on  the  power  of  the  X-ray  gen¬ 
erator. 

Observations  made  by  topographic  methods 
fall  into  three  groups—^locations  grown  into 
die  crystal,  dislocations  produced  during  plastic 
deformation  and  dislocations  introduced  or  modi- 
fiedby  division  processes.  The  lastof  these groups 
was  considered  by  in  Volume  1,  Point  Defects  (Pe¬ 
trenko  and  Whitworth  [1994]  see  section  onMolec- 
ular  Effects). 

Grown-in  dislocations 

All  normally  produced  crystals  (produced  with¬ 
out  special  precautions  b^g  taken)  or  grains 
widtin  polycrystals  contain  stable  networks  of 
grown-in  dislocations,  which  often  form  arrays 
OHistituting  low  angle  boundaries.  In  ice  almost 
allsudi  dislocationshave  diebasal  Burgers  vector, 
and  their  amcentrations  are  appreciably  lower 
than  those  found  in  metals  and  odier  commonly 
occurring  solids.  Examples  of  good  quality  ice 
from  die  Nfendoihall  Glacier  are  described  by 
Pukuda  and  Higashi  (1%9),  but  other  naturally 
occurring  ice  can  be  mudi  less  perfect  (Fukuda 
and  Sioji  1988).  For  die  study  of  individual  dislo¬ 
cations  try  topographic  methods,  mudi  lower  dis¬ 
location  densities  are  needed,  and  much  attention 
has  been  given  in  Japan  to  dw  growth  and  exam- 
inaticm  of  such  crystals  (Higashi  et  al.  1968,  Hi- 
ga^  1974,  Oguro  1988).  Single  crystals  diat  are 
several  centimeters  in  diameter  and  typically  10 
cm  Icmg  are  grown  in  glass  tubes  diat  are  seeded  by 
growdt  througli  a  capillary.  The  latest  refinement 
of  tills  method  is  that  of  CXitomo  et  al.  (1987).  Good 
crystals  contain  less  than  100  dislocations  with 
(a/3)  <2110>  Burgers  vectors  per  square  centime¬ 
ter,  but  there  are  usually  also  a  few  circular  or 
spiral  loops  that  have  [0001]  Burgers  vectors  and 
lie  on  or  close  to  the  basal  plane  (Fig.  7).  These 
loops,  which  are  of  prismatic  character,  have  been 
studied  by  Oguro  and  Higashi  (1981 ),  and  it  is  now 
believed  that  they  are  formed  by  the  condensation 
of  interstitials  (Oguro  et  al.  1988).  Crystals  with 
significant  concentrations  of  impurities  are  usual¬ 
ly  much  less  perfect  (Oguro  19^),  but  the  grains 
within  polycrystals  grown  slowly  under  carefully 
controlled  conditions  can  be  remarkably  good  (Liu 
etal.1992). 

Dislocations  with  plastic  deformation 

Mmyexperimentshavebeen  described  in  which 
dislocations  were  observed  to  move  and  to  multi- 


Figure  7.  Topographprojectedon  the  (0001) 
plane  showing  concentric  dislocation  loops 
inanas-grownaystalofice.  Thevariationof 
contrast  round  the  loop  is  an  in''vQtum  of  its 
prismatic  character.  The  othd  slocaHons 
seen  are  typical  of  the  low-density  random 
network  in  a  gpod  as-grown  crystal  (from 
Oguro  1988,  used  with  permission  cf  Hok¬ 
kaido  University  Press). 

ply  under  an  applied  stress.  Such  motion  repre¬ 
sents  an  extremely  early  stage  of  deformation; 
once  any  significant  macroscopic  strain  is  pro¬ 
duced,  foe  dislocation  density  becomes  too  high 
for  topographic  observations.  In  early  eiqieriments 
(e.g.,  Fukuda  and  Higashi  1973,  Jones  and  Gilra 
1973,  Mai  1976,  Fukuda  et  al.  1987,  Fukuda  and 
Higafoi  1988)much  information  was  lost  because, 
as  was  foimd  subsequently,  the  dislocation  struc¬ 
ture  changed  after  unloading  in  foe  time  foat  was 
required  to  obtain  foe  topographs.  In  more  recent 
work,  especially  that  using  synchrotron  radiation, 
it  is  much  easier  to  identify  di^rent  kinds  of 
dislocations  in  foe  topographs. 

Figure  8  is  a  topograph  showing  dislocations 
projected  on  foe  basal  plane  in  a  crystal  foat  had 
been  subjected  to  a  compressive  stress  in  a  vertical 
direction  (Ahmad  and  Whitworth  1988).  Features 
with  120°  angles  as  at  A  and  B  are  dislocations 
gliding  on  foebasal  plane.  Thesedislocations  glide 
as  almost  straight  segments  in  foe  screw  and  60° 
orientations,butfoecomersrapidlybecome  curved 
after  foe  stress  is  removed.  Very  occasionally,  for 
example  in  a  collapsing  loop,  basal  dislocations  in 
edge  orientation  have  also  been  seen.  The  long 
narrow  loops  like  those  at  C  and  D  have  basal 
Burgers  vectors  parallel  to  their  lengths  and  lie  on 
nonbasal  planes  oblique  to  foe  plane  of  t  rie  i:opo- 
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Figure  8.  Topograph  projected  on  the 
(0001)  plane  showing  dislocations  in¬ 
troduced  by  a  compressive  stress  pro¬ 
ducing  a  dtair  stress  on  thisplane  in  the 
verticttl  direction.  The  diffraction  vector 
is  1100;  A  and  D  are  lo^  in  the  basal 
plane.  CandDshowedgedislocation  seg- 
mentsonnonbasalplaneslyingoblufueto 
the  plane  of  the  topograph,  which  are 
dragging  long  screw  dislocations  bdiind 
them.  E  is  a  prismatic  loop  of  the  kind 
shown  in  Figure  7  (from  Ahmad  and 
Whitworth  1988). 


graph.  The  tip  of  the  loop  is  an  edge  dislocation  on 
the  ncmbasal  plane,  and  the  long  segments  in  the 
figure  are  screw  dislocations.  The  loop  E  is  a  pris¬ 
matic  loop  of  the  type  shown  in  Figure  7. 

Similar  features  to  the  nonbasal  loops  at  C  and 


D  have  been  reported  in  the  work  of  Fukuda  et  al. 
(1987)  (see  also  Fukuda  and  Higashi  1988),  but  the 
nature  of  this  nonbasal  glide  is  best  revealed  in 
crystals  stressed  parallel  to  the  basal  plane,  as 
shown  in  Figure  9  (^learwood  and  \^^twordl 
1989).  This  shows  dislocaticmspropagatingfroma 
scratch  on  the  back  surfece  of  a  crystal  with  die 
orientation  shown;  the  edge  dislocations  move  cm 
{ 1010}  planes,  but  the  screw  dislocaticms  diat  they 
leave  behind  do  not  glide  oa  diese  planes.  The 
separation  of  the  screw  segments  of  loops,  such  as 


Figure  9.  Sequence  of  topographs  (Uft)  sfwwing  edge  disloca¬ 
tions  propagating  on  the  (0110)  and  ( 1010)  planes  from  a 
scrat^  on  the  surface  ofa  crystal  oriented  as  Micated  above 
to  prevent  basal  slip  (from  Shearwood  and  Whitworth  1989, 
reproduced  courtesy  cfthe  Interrutional  Claciological  Soci¬ 
ety  from  Joimial  of  Glaciology,  35(120):  202). 


7 


C  in  Figiire  8,  arises  from  glide  on  the  basal  plane. 
Hoiulcd\  et  al.  (1990)  have  also  studied  nonbasal 
slip  in  specimens  with  the  [0001]  axis  perpendicu¬ 
lar  to  dve  stress. 

We  conclude  from  these  observations  that  all 
glide  dislocations  have  Burgers  vectors  in  the  bas- 
id  plane,  and  that  loops  expanding  cm  this  plane 
take  up  a  hexagcmal  form  made  up  of  screw  and 
60°  segments.  The  screw  segments  caimot  cross¬ 
glide  OR  n<mbasal  planes,  but  edge  dislocations 
can  glide  easily  on  nonbasal  planes  containing 
their  Burgers  vector.  Theseduuacteristks  are  prob¬ 
ably  unique  to  ke,  and  are  relevant  to  its  high 
degree  of  plastic  arusotropy. 


DISLOCATION  MOBILITY 

Experimental  observations 

The  most  detailed  measurements  of  dislocation 
velocities  are  those  of  Shearwood  and  Whitworth 
(1991),  who  obtained  sequences  of  topographs 
showing  the  positions  of  dislocations  between 
successive  applications  of  stress.  Figure  10  shows 
tracings  of  su^  positions  for  a  loop  expanding  on 
the  basal  plane,  and  Figure  1 1  shows  the  projectirm 
on  the  basal  plane  of  an  edge  disloca  ti<m  moving  to 
the  lower  left  on  a  nonbasal  plane.  For  the  stresses 
up  to  1  MPa  used  in  these  experiments,  the  dislo¬ 
cation  velocity  ua  was  found  to  be  directly  propor¬ 
tional  to  stress  t;  foe  velocities  per  unit  stress  for 
basalscrew,basal60°aiulnonbasaledgedisloca- 
ti(»is  in  pure  ice  are  plotted  as  functions  of  temper- 
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a.  Tracings  (f  the  positions  of  the  loop.  b.TopographicmagetracedasthesixthUnemFigurelOa. 

FigurelO.  Hexagonaldisloctttionloopexpandmgonthe{OOOl)plttnei^successiaeloadings(fromShettrwoodand 
Whitworth  1991). 


Figure  11.  Two  topographs  and  a  sequence  cf  tracings 
showingthenwtion^apinnteddislocationloopgtiding 
to  lower  left  on  a  nonbasal  plane  (from  Shearwood  and 
Whitworth  1991). 
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atureinHgutelZllieiesultsfornonbasaldisloca- 
tions  are  consistent  with  the  observations  of  Hon- 
doh  et  aL  (1990).  Eariier  results,  such  as  diose  of 
Yamamoto  and  Fukuda  (quoted  by  Fukuda  et  al. 
1987)  and  Mai  (1976),  do  not  distinguish  amcmg 
dideienttypes^dislocations,  and  were  probably 
subject  to  recovery  in  die  time  taken  to  obtain  the 
topographs.  Ihere  is  no  evidence  to  support  the 
nonlinear  dependenoeof  velocity  on  stressreport- 
edbyMaL 

In  so  far  as  die  Arrhenius  plots  in  Hgure  12  are 
strai^lines,dieacdvationeneigiesfordislocati(m 
glide  (actually  from  plots  of  ln(OdT/T)  vs.  1/7)  are 
as  in  Table  1. 

Anodier  technique,  which  in  principle  provides 
iriformationaboutdidocationmobility,is  internal 
diction.  Many  eqieiiments  (e.g.,  Vassoille  et  al. 
1978,TatibouetetaL  1986)haveob6ervedaoontri- 
butkxithat  increases  widiteirq>erature,isenhanoed 
by  plastic  deformation  and  often  depends  on  am- 


Figure  12.  Dislocatumvelocitiesper  unit  stress  vulxin 
pure  ice  as  /uruHons  of  inverse  tenqxrature  1/T  for 
screw  and  60“ dislocation  segments  on  the  basal  plane 
andforedgesegmentsonnonlHKidplttnes  (after  Shear- 
wood  and  WhUworA  1991). 

plitude.  This  edect  can  reliably  be  attributed  to 
dislocations  and  might  be  usable  for  distinguish¬ 
ing  amcmg  the  motion  of  kinks,  die  nudeation  of 
kinks  and  breakaway  from  pinning  points,  but 
there  are  toomany  disposable  parameters  for  us  to 
draw  condusions  from  die  oteervations. 

Peierls  model  for  basal  dislocations 
The  fact  that,  under  stress,  dislcxsitions  on  die 
basal  plane  glide  ({uickly  into  hexagcmal  form  and 
thoi  ^de  as  straight  segments  is  strong  evidence 
for  modem  across  a  Peierls  barrier,  as  is  also  ob¬ 
served  in  materials  such  as  silieem  (George  and 
Rabier 1987,  Nadgomyi 1988).  The  energy  of  a  dis- 
Icxaticmwillbeminiiriizedwhenitliesalongapar- 
ticular  line  in  die  crystal  lattice,  and  a  step  on  die 
dislcxaticm  between  two  such  lines  constitutes  a 
kink,  as  in  Figures  4  and  5.  The  kinks  can  glide 
alcmg  the  dislexation,  and  under  stress  their  mo¬ 
tion  will  carry  the  dislocaticm  forward  until  it 
ccHisists  of  almost  straight  segments  along  die 
directions  of  minimum  energy.  In  ice  diese  are  the 
screw  and  60°  segments  seen  in  Hgure  10.  To  ad¬ 
vance  fiirdier,  double  kinks  must  be  thrown  for¬ 
ward  across  die  Peierls  barrier  and  then  move 
sideways  until  they  reach  die  end  of  die  segmait 
or  are  annihilated  by  other  kinks.  A  computer 
simulaticm  illustrating  diis  bdiavior  is  shown  in 
Hgure  13.  The  dislcKaticm  velocity  is  given  by 

Pd  =  foikOk  (1) 

where  Hk  .  number  of  kinks  per  unit  langdi 


F^rel3.  Computer  shnulatumcf glide  cfadisbcation  across  a  Peierls  barrier 
bp  the  nudeatkm  cfkmkpairs  and  by  the  glide  tfthe  kinks  along  the  dislocation. 
The  dislocation  is  gliding  upwards  in  the  figure  and  successive  positions  are 
shown  dbooe  one  aruither.  The  original  positkmfbr  each  case  is  indicated  by  the 


marks  at  the  edges  cfthe  diagram. 

=  their  velocity 

h  =  separatitm  of  ttie  Peieris  trough. 

According  to  tfie  standard  theory  of  tfiis  kind  of 
diskxatkxt  motion  (Hirtti  and  Lo^  19S2),  in  die 
liinitofk>w8tre88(iM«ji;bT,whereaisthelattice 
parameter  andbis  die  Burgers  vector),Od  is  linear- 
ty  prc^xntional  to  stress.  It  depends  on  bodi  the 
activation  free  enei:^  to  form  an  isolated  kink  Fk 
and  duit  to  move  it  Fq,  according  to  die  equaticai 

oj.vZlh^exp .  (2) 

The  quantity  v  is  a  diaracteristic  frequency  that 
caimot  be  more  than  the  Ddiye  cutoff  freqii^cy. 
The  free  energy  Fv  can  be  written  in  terms  of  the 
energies  and  entropies  of  formation  [£|((T)  -  75|J, 
and  diemeasuied  activation  energy  for  glide  is  the 
teiiqperatureindependentpart<d(E|(+  £„).  Howev¬ 
er,  widi  the  above  eiqierimental  values  and  other 
known  parameters,  tteentn^nesof  activation  turn 
out  to  lx  remarkabty  high  (for  details  see  9iear- 
woodand  Whitwordi  1991).  This  model  indicates 
diat  dislocations  seem  able  to  move  more  easily 


than  expected.  The  contributions  arising  from  kink 
nudeation  and  Idnk  migration  can  m  principle  be 
separated  by  studying  die  motion  of  curved  seg¬ 
ments  (Hoiididi  1992),  but  diis  has  not  yet  been 
successfully  achieved. 

Protim  disorder 

Aunk]uefeatureofice,  first  recognized  by  Gkn 
(1968),  is  diat  the  discvder  of  the  protcms  presents 
in  principle  an  obstacle  to  die  glide  of  dislcxadcms. 
This  arises  quite  simply,  because,  if  two  planes  of 
molecules  are  linked  by  randomly  orientedhydro- 
genbonds,  they  carmotbe  sheared  over  one  anoth¬ 
er  and  still  link  up  correcdy.  The  idea  is  easily 
understood  by  considering  a  60°  dislocaticm  on 
planes  of  the  shuffle  set  as  shown  in  Hgure  14, 
where  die  positions  of  die  protons  are  disordered. 
For  the  dislcxation  to  glide  to  the  left,  bond  DIX 
must  be  broken  and  U  joined  to  C.  This  presents 
no  problem,  but  for  motion  to  die  right,  C  must 
link  to  C,  which  would  create  a  D-defect.  There  is 
not  sufficient  energy  from  the  stress  to  create  this 
defect,  and  there  is  in  general  no  local  rearrange¬ 
ment  of  bonds  that  will  avoid  die  frnmation  of  a 
defect  somewhere.  In  practice  we  expect  glide  to 


Figure  14.  Section  in  the  (1 100)  plane  cfa  60” dislocation  on  a 
plane  of  the  shuffle  set  in  the  structure  of  ke,  Qlustrating  the 
consequence  of  the  disorder  cf  the  protons  according  to  the  modd 
of  Glen  (1968). 


occur  by  demotion  of  kinks  along  the  dislocation, 
widi  eadi  step  of  die  kink  involving  a  single  ex¬ 
change  ofbondslike  the  one  just  oxisideted.  There 
will  be  a  50%  diance  ai  the  bonds  being  mis- 
matdied  at  eadi  step. 

Glen  therefore  prc^iosed  foat  the  rate-limiting 
step  for  dislocation  motion  may  be  ttie  rate  at 
tdiidi  bonds  are  randomly  reoriented  by  ions  or 
Bjerrum  defects,  as  m  the  process  of  dielectric 
relaxation.  It  is  important  to  realize  duit  die  stress 
cannot  foroediereocientation  of  therequiredbonds 
asthedislocBtkmapproaches.'niisideawasquan- 
tified  fm  kinks  on  die  shuffle  plane  by  Whitwordi 
et  aL  (1976)  and  Frost  et  aL  (1976).  The  stresses 
invtdved  are  always  such  diat  diey  impose  a  small 
bias  on  dw  random  motion  of  die  kitila;  die  kiiiks 
arenot  pushed  up  against  dw  mismatched  bonds. 
Widi  dlls  assumption,  die  kink  velocity  for 
reorientations  by  Bjerrum  defects  is  given  by 


Tx^fisT 


(3) 


where  is  die  mean  time  to  reorient  a  bond. 

For  dislocations  on  planes  of  die  glide  set,  the 
obstacle  presented  by  proton  disorder  may  be  less 
severe.  Ptovided  die  cxne  is  not  reconstructed, 
molecules  B  and  C  in  Hgures  5a  and  b  are  free  to 
rotate  about  die  sin^  bond,  perpendicular  to  die 
glide  plane,  as  diey  switch  dieir  linkages  from  G 
and  H  to  A  and  D.  For  a  kink  on  an  unreconstruct¬ 
ed  partial  dislocation,  Whitwordi  (1980)  showed 
diat 


p.  glZffiVSil*  t.  (4) 

If  die  partial  dislocatkm  is  reconstructed,  the  bar¬ 
rier  fnesented  by  proton  dlscmier  will  be  mudi 


greater,  but  this  situation  has  not  been  analyzed 
dieoretically.  In  all  cases  where  a  Peierls  barrier  is 
present,  proton  disorder  will  affect  the  rate  kink 
nudeation  as  well  as  die  Idnk  mobility,  but  it  can 
only  make  die  theoretical  velocity  lew  dian  that 
predicted  from  eq  2.  A  dieory  has  also  been  devel¬ 
oped  for  a  flexible  dislocation  line  (Whitworth 
1983).  Ihis  does  not  seem  to  be  applicable  to  dislo¬ 
cations  that  glide  as  straiglit  segments,  but  may  be 
applicable  f^  nonbasal  edge  dislocations. 

Ihe  quantity  q,  is  die  mean  time  for  the  reorien¬ 
tation  of  bonds  dose  to  die  dislocation  core,  and,  if 
we  assume  diat  die  ice  rules  axe  applicable,  diis 
will  take  place  by  the  motion  of  Bjerrum  defects  or 
ions.  Inb^crystal,  diis  time  isappxoximately  die 
same  as  die  Dd>ye  relaxation  time  but  Shear- 
wood  and  Whitwordi  (1991)  have  shown  diat  for 
dislocations  to  move  at  die  observed  rate,  die 
appropriate  value  of  q,  must  be  much  shorter  dian 
diik  A  critical  experiment  to  establish  the  possible 
relation  of  the  dislocation  velodty  to  die  bond 
reorientation  rate  in  bulk  ke  is  to  measure  velod- 
ties  in  doped  ke.  Hus  was  first  attempted  for  HF 
dealing  by  Mai  et  aL  (1978),  who  found  a  small 
efi^,  but  mudi  less  dum  predicted.  Using  HQ 
doping,  Sheaxwood  and  Wliitwordi  (1992)  found 
no  significant  effect  under  conditions  where  qj 
was  known  to  have  been  reduced  by  more  dian  a 
factor  of  10. 

From  all  diat  we  know  about  the  disorder  of  the 
protcHis  in  ke,  it  is  essential  diat  diere  be  a  process 
for  die  reorientatirm  of  bonds  at  the  dudocation 
core  duit  retains  compatibility  with  the  ice  rules  in 
die  surrounding  material.  It  could  arise  frcrni  an 
enhanced  concentration  or  mobility  of  icais  or 
Bjerrum  defects  near  the  core.  Interstitials  or  va¬ 
cancies  should  not  be  directly  involved  as  they  do 
not  change  die  proton  disorder.  Perez  et  aL  (1978, 
1980)  postulated  diat  die  dislocatkm  core  was  in 
some  sense  noncrystalline,  thereby  avoiding  die 
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obstacte  presented  by  ttie  ke  rules.  Their  particular 
model  has  several  adjustable  parameters,  and  these 
take  unphyskal  values  diat  depend  on  the  nonlin¬ 
earity  in  dieir  observations  of  The  fact  that 
dislocatkms  glide  as  straight  segments  in  crystal¬ 
line  orientations  indicates  that  the  core  must  retain 
mudi  of  its  crystalline  character.  The  precise  way  in 
which  dislocations  overcome  the  barrier  presented 
by  proton  disorder  is  at  present  unknown,  and,  as 
wewiUseemdte  Doped  Ciystttlsand  Electrical  Effixts 
section,  there  is  conflicting  evidence  about  whedter 
proton  disorder  is  rate  limiting. 

There  is  no  evidence  for  a  Peierls  barrier  to  the 
motion  of  ruxnbasal  edge  dislocations,  and  this  is 
consistent  widi  dieir  haviirg  a  lower  activation  en¬ 
ergy  fw  glide.  However,  die  limitations  presented 
by  proton  disorder  are  in  principle  equally  impor¬ 
tant  in  diis  case  (Whitworth  1983,  Shearwood  and 
Whitworth  1991).  The  most  significant  flK^t  about 
nonbasal  slip  is  the  complete  absence  of  any  glide  of 
screw  dislocations  on  the  nonbasal  platin  at  die 
stresses  used.  This  is  strong  evidence  duit  the  screw 
dislocations  are  dissociated  into  partial  disloca¬ 
tions  on  die  basal  plane,  and  provides  circumstan¬ 
tial  evidence  for  basal  slqi  beteg  between  planes  of 
die  glide  set 

Weertman  (1963)  proposed  a  furdier  mechanism 
by  whkh  die  disoirder  of  the  protons  in  ice  will 
inhibit  die  nukian  of  dislocations:  anelastk  loss, 
iriikh  arises  from  reorientation  of  molecules  in  die 
stress  field  of  die  dislocation  as  it  moves.  Any  such 
efiect  will  constitute  a  barrier  additional  to  those 
aheadyconsidered,  and  seems  likely  to  be  omipara- 
tivdy  smalL 

ROLE  OF  DISLOCATIONS  IN 

plastic  deformation  of 

SINGLE  CRYSTALS 
Pore  ct]fata]s 

For  die  plastic  deformation  of  a  crystal  on  a 
sing^  sl4>  system,  die  strain  rate  e  is  given  by  the 
equation 

e«#NdliPd  (5) 

whereNd  «  density  of  mobile  dislocations  (length 
per  unit  vcdume) 
b  3  Burgers  vector 
Ud  *  dislocation  velocity 
^  «  foctor  depending  cm  the  orientation  of 
die  sl4>  system. 

If8everalsltysystemsoperatetogedier,dieirc(mtri- 
butions  to  strain  rate  must  be  summed.  Crystals 


of  ice  normally  deform  almost  exclusively  on  the 
basal  plane,  and  Figure  15  shows  the  form  of  the 
stress-strain  curve  for  such  deformaticm  in  ccm- 
stant  strain  rate  taisile  tests  (Higashi  et  aL  1964). 
The  peak  followed  by  a  yield  drop  arises  because 
Nd  is  initially  small  and  the  ccmstant  deformation 
rate  requires  a  high  value  of  but  as  increases 
during  the  deformation,  the  stress  necessary  to 
maintain  a  smaller  t>d  becomes  smaller.  Evoitual- 
ly,  deformation  may  pnxeed  at  constant  stress. 
There  is  little  work-h^ening  in  ice,  and  in  die 
steady-state  region  there  will  be  a  balance  among 
dislocation  multiplicaticm,  emergence  from  die 
surfoce  and  recovery  processes. 


Figure  15.  Stress-strain  curves  for  tensile 
d^ormatUmofsingle  crystals  of  ice  at -15X^ 
and  at  the  constant  strain  rates  shown  (c^ter 
Higashi  et  dl.  1964). 

The  correspcmding  behavior  to  Figure  15  in  a 
creep  test  at  constant  stress  would  be  an  initially 
accelerating  creep  leading  to  a  constant  strain  rate. 
Higashi  et  al.  (1965)  have  reported  this  kind  of 
creep  in  a  very  special  bending  geometry,  but  in 
tensile  tests  ccmtinucnisly  accelerating  aeep  is 
observed,  for  which  e  ^  widi  m  ■■  15  to  2.  The 
macroscopic  deformaticm  of  ice  depends  critically 
cm  die  initial  state  of  die  crystal  anci  die  ccmditicms 
of  the  experiment,  but  results  are  gaierally  fitted 
to  die  empirical  equaticm 

e  «*  o’*  exp  (-  E/koT)  (6) 

where  for  single  crystals  n  »  2.  There  is  some  dif¬ 
ficulty  in  determining  a  value  for  die  activation 
energy  £,  because  values  deduced  fiom  ccmstant- 
strain-rate  tests  depend  cm  the  value  of  n,  while 
strictly  comparable  ccmditicms  are  difficult  to 
achieve  in  creep  tests  at  difierent  temperatures. 
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Td>lc  2.  AclivaticMi  cnogics  for  plastic  deformation 
of  sin^c  crystals  of  ice  at  temperatures  of  -10  to 

-«rc 


iehnnttion 

Actamtian 

emrgjfE 

(eV) 

Reference 

Tensik-«anstant  f 

0.6S 

Higashi  et  aL  (1964) 

Teraile-caratant  E 

0.62 

Readey  and  Kingety  (1964) 

Bending  deep 

0.68 

Higashi  et  aL  (1965) 

Temiie  creep 

0.68 

Jones  and  Glen  (19^) 

Tenaile  creep 

aso 

Homer  and  Glen  (1978) 

Tensile  creep 

0.62 

Ramaeier  (unpublished)* 

*  Fran  Weertman  (1973). 

Table  2  summarizes  values  deduced  by  different 
authoisindifferentwaysfromtestsonsinglecrys- 
talsindierange-lOto-SO^C.TlMrearemanymore 
data  on  polyoystalline  ice^  but  diis  deforms  very 
difierently  and  other  omsiderations  will  be  rate 
limiting,  ^r  a  detailed  review  of  creep  in  ice,  see 
Weertman  (1973). 

The  interesting  feet  that  emerges  from  this  table 
is  that  the  values  of  £are  all  less  than  those  in  Table 
1  for  foe  glide  of  dislocaticm  segments  on  die  basal 
plane  in  diis  range  of  temperatures.  Extensive 
studies  of  silicon  reveal  no  sudi  discrepancy,  but 
it  does  exist  in  some  ofoer  semiconductors  (Rabier 
and  George  1987).  There  are  two  possible  reasons 
for  die  difierence.  Eidier  foe  dislocation  density 
for  steady-state  creep  increases  with  falling  tem¬ 
perature,  or,  under  deformation  conditions,  foe 
dislocatifms  do  not  glide  as  straight  segments.  If  it 
is  easier  to  create  kinks  widiin  foe  deforming 
material  diantmder  the  near-perfect  cemditions  of 
topography  experiments,  the  activation  energy 
shxild  be  reduced.  Macroscopic  deformation  is 
sudi  a  complicated  process  diat  it  is  not  surprising 
diat  diere  is  ccmsiderable  variation  among  exper¬ 
iments. 

Crystals  oriented  wifo  die  stress  perpendia^ 
to  die  [(XX)1]  axiscanbe  deformedby  slip  on  ( IlOO) 


planes,  but  this  requires  a  stress  of  some  50  times 
that  for  slip  on  the  basal  plane,  and  after  such 
deformation  voids  are  created  in  the  specimen 
(Muguruma  et  al.  1966).  We  have  seen  in  topo¬ 
graphic  experiments  diat  edge  dislocations  move 
very  easily  on  these  nonbasal  planes,  but  screw 
dislocations  were  not  observed  to  glide  off  the 
basal  plane.  This  means  that  a  dislocation  loop  that 
starts  to  move  on  a  ncmbasal  plane  cannot  expand 
to  cover  foe  whole  plane;  slip  will  be  confined  to 
narrow  strips  bounded  by  screw  dislocations  as 
seen  in  Hgure  9,  and  macroscopic  slip  will  be  very 
difficult.  Hondoh  et  al.  (1990)  reproduce  a  figure 
showing  that  such  slip  is  associated  with  very 
foort  slip  lines,  which  correqxnid  well  with  this 
interpretation. 

To  produce  macroscopic  deformation,  it  is  not 
sufficient  that  dislocatimis  should  glide.  There 
have  to  be  processes  by  which  dislocations  can 
multiply  on  tiieir  slip  plane  and  by  whidi  slip  can 
be  transferred  from  one  slip  plane  to  anotiier.  The 
standard  mechanism  is  ihe  Frank-Read  source 
(see  Read  1953).  An  example  of  such  a  source  is 
illustrated  in  Figure  16a,  in  which  dislocation  seg¬ 
ments  spiral  around  a  fixed  point  where  foe  dislo¬ 
cation  makes  a  step  from  one  glide  plane  to  anoth¬ 
er.  Ahmad  et  al.  (1986)  observed  sudi  a  source  in 
topographic  experiments  on  ice,  and  Figure  17 
foows  an  example.  In  normal  materials  many  sudi 
sources  are  generated  by  the  cross  slip  of  screw 
dislocations  off  tiieir  primary  glide  planes,  but  tiie 
immobility  of  screw  dislocations  on  foe  nonbasal 
planes  in  ice  means  that  this  does  not  occur.  How¬ 
ever,  a  feature  of  ke  is  foe  high  mobility  of  edge 
segments  on  nonbasal  planes,  and  this  means  tiiat 
foe  segment  S,  whidi  acts  as  foe  fixed  center  of  foe 
Frank-Read  action  in  Hguie  16a,  may  not  remain 
fixed  but  more  often  glides  away  from  the  original 
dislocation,  trailing  screw  segments  behind  it  as 
shown  in  Figure  16b.  This  leads  to  features  sudi  as 


FigurelS.  OpenOkmi^aFnmk-Readsourceon  thebasdplane(a),(mdvjhathappens 
inkx^thetinkmgs^mentSdoesnotrmttmfixedbuiglidestttpuilyonthemndMisal 
slippUme(b). 
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Ftgurel?.  Setiuencecftopographspro]ectedonthe{(XXn)planeshomngoperattonofttFrank-Readsourcemice.  This 
corresponds  to  Figure  16a  (from  Aknud  et  al.  1992,  used  with  permission  of  Hokkaido  Unixxrsity  Press). 
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Fipae  18.  Sequence  of  topographs  projected  on  the  (0001)  plane  showing  dislocation  multiplkation 
arisingfrom  diefastedgesegmentssuchasAonanonbasalplaneobliqueto  theplaneofthefigure.  The 
sectkmidB,wt^devdopshexttgpnalfeatures,ison  the  ba^plane  and  glides  more  slowbf.  In  thelast 
two  topographs  the  dislocation  loops  have  cut  the  surface  at  upper  left  (from  Ahmad  et  al.  1992,  used 
withpermission  tfHoUcakb  Unioersity  Press). 


Tune  (h)  - ► 

a.  Creep  curves. 


Resolved  strain  (%)  — ► 

b.  Constant-strain-rate  tests. 


Figure  19.  Test  results  from  pure  single  crystals  cfitx  and 
crystals  doped  with  HF  at  the  concentrations  gmen.  Temper¬ 
ature  is  -70°C  (after  Jones  1967). 


C  in  Hgure  8  and  to  multiplication  of  the  kind 
shown  in  Figure  18.  The  behavior  of  nonbasal  edge 
dislocations  that  cut  die  surface  is  also  very  signif¬ 
icant  (Ahmad  et  al.  1992,  Shearwood  and  Whit¬ 
worth  1993). 

Doped  crystals  and  electrical  effects 

Jones  and  Glen  (1969b)  found  that  doping  ice 
widi  HF  produced  a  remarkable  softening  in  both 
creep  tests  and  constant-strain-rate  tests  at  -60  to 
-70®C(Fig.  19).  The  softeiing  was  observed  both  in 
crystals  doped  before  deformation  and  when  HF 
was  diffused  into  the  specimen  part  way  through 
die  test  NH40Hproducedasmallhardening.Naka- 
mura  and  Jones  (1970)  deformed  HCl-doped  ice, 
but  at  high»  temperatures,  and  observed  a  soften¬ 
ing,  diough  the  edect  was  much  less  than  with  HF. 
All  of  these  observadcms  indicate  a  strong  corre- 
latum  in  the  effects  of  die  dopant  on  die  deforma¬ 
tion  rate  and  cm  the  rate  of  bond  reorientation  as 


observed  in  dielectric  experiments.  This  led  to 
Glen's  (1968)  hypothesis  that  proton  disorder  is 
the  rate  limiting  process  for  dislocation  glide. 

The  topographic  observations  of  Shearwood 
and  Whitworth  (1992)  showed  no  effect  of  HCl 
doping  on  the  mobility  of  dislocations  at  a  level  of 
doping  that  produced  a  ten-fold  decrease  in  To, 
but  this  doping  was  very  much  less  than  that 
which  gave  the  large  effects  observed  by  Jones  and 
Glen.  A  further  complication  is  that  strain  rate  de¬ 
pends  cm  the  product  of  dislocation  density  and 
dislocation  velocity,  and  Jones  and  Gilra  (1972) 
found  that  diffusing  HF  into  ice,  as  was  done  in 
some  of  the  experiments  referred  to,  produced  a 
large  increase  in  dislocation  density. 

An  alternative  way  of  changing  the  concentra¬ 
tions  of  point  defects  that  gives  rise  to  bond  reori¬ 
entation  is  the  application  of  an  electric  field  (Pe¬ 
trenko  and  Schulson  1992).  By  applying  electric 
fields  to  tiiin  specimens  undergoing  deformation 
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in  shear,  Petrenko  and  Schulson  (1993)  have  shown 
that  a  reduction  in  the  high-frequency  conductivity 
leads  to  a  corresponding  reduction  in  creep  rate. 
We  have  already  seen  that  dislocation  glide  in¬ 
volves  the  reorientation  of  bonds  near  the  core. 
There  is  coirsiderable  evidence  that  under  certain 
conditions  changing  the  bond  reorientation  rate  in 
the  bulk  of  the  crystal  can  have  an  effect  on  the 
dislocation  mobility. 

Different  kinds  of  electrical  effects  can  arise  if  the 
dislocations  carry  a  net  charge.  Such  a  charge  is  to 
be  expected  where  there  are  dangling  bonds,  be¬ 
cause  the  numbers  of  bonds  with  and  without 
protons  do  not  have  to  be  equal.  In  equilibrium  any 
such  charge  will  be  screened  by  a  siurounding 
cloud  of  excess  electrical  point  defects  of  opposite 
sign  (see  Whitword)  1975),  but  during  deformation 
adislocationmaybecome  separated  from  discharge 
cloud.  Petrenko  and  Whitworth  (1983)  observed 
small  transverse  electric  currents  associated  with 
the  tensile  defbrmaticm  of  previously  bait  crystals, 
and  interpreted  them  as  asoibable  to  dislocations 
carrying  a  net  positive  charge  of  at  least  0.002  pro¬ 
tonic  charges  per  atomic  length.  Itagald  (1970)  has 
reported  X-ray  topographic  experiments  in  which 
dikocations  appear  to  move  in  an  alternating  elec¬ 
tric  field  In  cases  where  the  sign  of  the  cnarge  could 
be  identified,  it  was  positive,  but  it  is  difficult  to 
deduoemagnitudesftomsuchexperiments.Aprob- 
lem  with  any  experiment  of  this  land  is  that  the  field 
should  be  maintained  for  a  sufficient  time  to  move 
die  dislocations  by  an  observable  amount,  but  in 
this  time  die  field  is  largely  eliminated  by  polariza¬ 
tion  and  ctxiduction  in  the  ice.  If  die  ^location 
cores  are  indeed  reomstructed,  as  suggested  in  the 
Basal  Dislocations  section,  then  any  core  charge  will 
be  confined  to  kink  sites  or  perhaps  to  ^mim  type 
defects  trapped  in  the  reconstruction. 

STACKING  FAULTS 

Structure  of  stacking  fiiults  in  ice 

The  crystal  structure  of  ice  consists  of  (0001) 
planes  of  molecules  stacked  above  one  anoffier  as 
illustrated  in  Hgure  2.  If  die  positions  of  diese 
layers  in  their  planes  are  denoted  by  the  letters 
A,B,C,  as  defined  in  Hgure  3,  die  hexagonal  struc¬ 
ture  of  ice  Ih  follows  die  sequence 

AABBAABBAABB.... 

TheCpositions  are  unoccupied,  resulting  in  empty 
diann^runningdiroug^^sttuctureindie[0001] 
dkectkm.  Stacking  fiiults  arise  if  die  stacking  across 


the  planes  marked  "glide"  in  Figure  2  departs 
from  the  perfect  sequence.  For  a  general  diknis- 
sion  of  such  faults,  foe  reader  is  referred  to  Hirth 
and  Lothe  (1982)  or  a  similar  text.  To  simplify  the 
discussion,  we  will  denote  each  pair  of  layers  such 
as  AA  by  the  single  letter  A,  and  the  stack-ing  is 
then  that  familiar  in  hexagonal  close  packing 

ABABABAB... . 

Cubic  ice  (ice  1^)  has  the  stacking 
ABCABCABC... 

but  is  energetically  unstable  relative  to  ice  Ih- 
A  stacking  fault  is  a  planar  defect  normally  ly¬ 
ing  on  a  (0001)  plane  of  the  glide  set.  It  must  extend 
to  the  surface  or  terminate  at  a  partial  dislocation 
widi  a  Burgers  vector  equal  to  die  displacement 
required  to  create  the  fault.  To  distinguish  among 
the  possible  kinds  of  faults,  it  is  useful  to  introduce 
Frank's  notation,  which  ciHicentrates  not  on  the 
absolute  location  of  the  layers  A,B,C  but  only  on 
the  stacking  relative  to  the  layer  below.  Ihus,  the 
equivalent  stackings  of  B  on  A,  C  on  B  and  A  on  C 
are  all  denoted  by  the  symbol  A,  while  the  inverse 
stacking  of  A  on  B,  B  on  C  or  C  on  A  are  given  the 
symbol  V  The  hexagonal  stackings  ABABAB..., 
BCBCBC...  or  CACACA...  are  then  all  denoted  by 

AVA7AVAV... 

whereas  cubic  stacking  would  be 

AAAAAAAA...  or  WWVWV... . 

There  are  four  simple  ways  of  introducing  stack¬ 
ing  fiiults  into  ke,  and  these  are  illustrated  in 
Figure  20.  The  first  is  to  shear  a  B  layer  over  an  A 
layer  into  a  C  position.  This  changes  a  A  into  a  V 
giving  the  fiiult 

i 

BABACBCB...  . 

VAWVAV 

This  fiiult  is  iUustrated  in  Figure  20a,  in  which  it 
terminates  in  a  partial  dislocation  of  the  Burgers 
vector  (af3)  <  1010>.  When  a  perfect  (afi)  <  1 120> 
dislocation  dissociates  into  two  9iockley  partial 
dislocaticHis  on  a  (0001)  plane  of  the  ^de  set 
according  to  the  equation 

(a/3)  <1120>  (a/3)  <10T0>  +  (a/3)  <0lT0> 


c.C  layer  added.  d,Clayeraddedfbttowedbyshearabovetheaddedlayer. 

Figure  20.  Projections  (tf  structure  of  ice  Ii,ona  (1210)  plane  ^wwing  stacking  faults  in  the  right-hand  portion  of 
each  dtagram  terminating  at  an  appropriate  dislocation  in  the  center. 


as  described  in  the  Baso/ Oiis/ocafions  section,  tfiis  is 
tiie  type  of  stacking  tiniltribbon  that  will  be  formed 
between  tiiem. 

Hie  other  types  of  faults  require  the  addition  or 
removal  of  alayer  A,  B  or  C,  which  ccmsists  of  two 
planes  of  molecules.  If  an  A  layer  is  removed,  the 
B  layers  on  opposite  rides  of  it  can  only  link  to- 
gettier  if  there  is  a  displacement  of  one-half  ofttie 
crystal  over  tiie  other  by  an  amount  (a /3)  <  1010>, 
generating  the  foult 

i 

ABABCBCB... 

VAVWAV 


This  fault  is  shown  in  Figure  2(%i.  Hie  dislcKation 
surrounding  it  has  both  prianatic  and  glid^com- 
ponents,  and  its  Burgers  vector  is  (1/6)  <2023>. 

Unlike  the  case  just  described,  the  addition  of  a 
C  layer  to  tiie  perfect  ABAB. . .  structure  does  not 
necessitate  any  shear,  and  the  fault  generated  is 

i 

ABABCABCB... 

AVAAAAVA 

Hiis  fault  is  illustrated  in  Figure  20c;  it  is  sur¬ 
rounded  by  a  prismatic  dislcKaticm  loop  of  Burg¬ 
ers  vector  (c/2)  [0001].  The  fault  contains  four 
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consecutive  A-type  stackings,  which  resembles 
cubkkeJtcan  lower  itsaiergybyashear  between 
the  C  and  die  A  layer  at  the  point  marked  i. 
yielding  die  fault 

i 

ABABCBCBC... 

AVAAVAVA 

which  is  shown  in  Figure  20d.  Examination  of  the 
AV  sequence  shows  that  this  is  the  same  stacking 
as  in  Hgure  20b,  but  the  nature  of  the  bounding 
dislocatimi  corresponds  to  an  interstitial  rather 
dianavacancy-typeprismaticloop.Stacking  faults 
are  often  dashed  as  tntrmstc  if  perfect  stacking  is 
maintaineduptotheplaneofthehiultandextnn- 
stc  if  it  is  not;  on  diis  basis  the  faults  of  Figures  20a, 
b  and  d  are  intrinsic  and  only  the  higher  energy 
hiult  20c  is  extrinsic. 

Observations  of  stacking  faults 

Stacking  faults  of  macroscopic  dimensions  can 
be  observed  by  X-ray  topography  and  have  been 


studied  extensively  in  ice  by  the  group  at  Sapporo. 
Their  work  is  reviewed  by  Fukuda  et  al.  (1987), 
Oguro  and  Hondoh  (1988)  and  Oguro  et  al.  (1988). 
Stacking  foults  in  ice  are  energetic^y  unstable  and 
never  remain  in  crystals  that  have  been  well  an¬ 
nealed  and  observed  at  the  temperature  of  anneal¬ 
ing.  However,  a  few  large-area  faults  do  occur  in 
fr^hly  grown  crystals  (Hondoh  et  al.  1983),  and 
their  formation  is  enhanced  by  doping,  particularly 
with  NH3  (Oguro  and  Higashi  1973). 

Figure  21  is  an  example  of  a  topograph  of  pure 
ice  cooled  to  -45^C.  The  dark  areas  correspond  to 
stacking  faults  parallel  to  die  basal  plane.  The  large, 
irregular  faults  originate  from  the  growth  of  die 
crystal,  while  the  smaller  patdies  are  prismatic 
loops  formed  during  cooling.  The  presence  of  con¬ 
trast  within  the  area  of  the  fault  shows  that  the  fault 
vector  f  satisfies  the  condition  g-f  ^  0  or  an  integer, 
where  g  is  thediffraction  vector  of  die  topograph. 
With  g  s  <1010>  this  means  duit  the  faults  must 
have  a  shear  compment.  The  bounding  disloca- 
timishave  a  [0001]  component,  and  careful  analysis 
shows  that  Ae  fiiulted  loops  formed  on  cooling  are 


FigureZl.  X-ray  topographpro}ectedma{0001) 
plane  of  a  crystal  of  pure  ice  cooled  rapidly  to 
-45°C,  resulting  in  the  production  of  fdidted 
prismatic  dislocation  loops  by  the  condensation 
ofinterstitials(fromOguroetal.  1988,usedmth 
permission  of  Hokkaido  University  Press). 


interstitial  in  character  (Hondoh  at  al.  1983).  It 
appears  that  an  interstitial  loop  with  the  Burgers 
vector  (c/2)  [0001]  does  in  fact  lower  its  energy  by 
taking  the  sheared  form  of  Figure  20d  rather  than 
that  of  Figure  20c.  The  stable  prismatic  loops  shown 
in  Figure  7  are  presumably  perf^  dislocations, 
with  the  Burgers  vector  c[0001]  including  two  dou¬ 
ble  layers  of  molecules  within  the  loop;  such  dislo¬ 
cations  will  have  high  energy  and  will  almost  cer¬ 
tainly  be  dissociated,  but  the  stacking  favdt  ribbon 
between  them  is  unresolvable  in  topography. 

As  described  in  Volume  1,  Point  De^s  (Petren¬ 
ko  and  Whitworth  [1994],  see  section  on  Molecular 
Defects),  very  detailed  studies  of  the  growth  and 
shrinkage  of  hiulted  and  imhiulted  dislocation  loops 
have  yielded  the  best  available  parameters  for  the 
formation  and  self-diffusion  of  interstitials  in  ice 
(Goto  et  el.  1986,  Hondoh  1992).  The  energy  of  a 
(1/6)  <  2023>  fault  has  been  determined  to  be  031 
mj  m“2,  which  is  3  x  10^  eV  per  molecule  in  the 
plane.  This  is  assumed  to  be  the  lowest  energy  type 
of  foult  because  it  has  oitly  two  adjacent  A's  in  the 
stacking  sequettce.  The  shear  fault  with  diree  A's  is 
estimated  to  have  twice  this  energy,  and  on  this 
basis  the  separation  of  the  9iockley  partial  disloca¬ 
tions  fbrm^  by  the  dissociation  of  a  perf^  (a/3) 
<1120>  dislocation  on  a  basal  plane  has  been  esti¬ 
mated  by  Hondoh  et  al.  (1983)  to  be  20  run  for  a 
screw  dislocaticm  and  46  run  for  an  edge  disloca¬ 
tion.  Observaticm  of  such  dissociation  lequireselec- 
tron  microscopy  and  has  not  yet  been  achieved  in 
ice. 

GRAIN  BOUNDARIES 
Structure 

Natural  ice  is  polycrystalline  and  theboundaries 
between  the  individual  crystals  (or  "grains")  canbe 
ccmsidered  as  planar  defects  within  otherwise  per¬ 
fect  material.  We  will  briefly  describe  the  intrinsic 
properties  of  ttiis  class  of  defects,  but  will  not  here 
amsider  foeir  role  in  ttie  nuicroscopic  properties  of 
polycrystallme  ice. 

Examinatian  of  polyaystalline  samples  shows 
grains  of  many  ^pes,  sizes  and  orientations,  de¬ 
pending  tm  the  history  of  the  material  (e.g.,Matsu- 
da  and  Wakahama  1978).  In  well-annealed  ioe,  the 
boundaries  are  fairly  i^t,  which  is  a  necessary 
conditiiHi  for  minimizing  foe  surfece  energy.  With¬ 
in  a  single  grain,  or  within  a  piece  of  ke  that  is 
nominallya8ingleciystal,foerewilloftenbesub- 
boundaries,acro6swhichfoereisalatticemisorien- 
tation  of  at  most  a  few  degrees  and  often  very  mudr 
less.  These  boundaries  are  made  up  of  arrays  of 


dislocations  in  accordance  with  geometrical  prin¬ 
ciples  that  apply  to  any  crystalline  material  (see  for 
example  Read  1953  or  Hirth  and  Lothe  1982). 
Except  for  elastic  interactions  with  one  another  to 
form  a  stable  structure,  these  dislocations  behave 
as  individual  dislocations  within  a  single  crystal. 
A  suitable  sub-boimdary  will  migrate  by  foe  glide 
of  these  dislocations  under  an  appropriately  ori¬ 
ented  stress,  as  has  been  observed  in  ice  by  Higashi 
and  Sakai  (1961). 

For  large-angle  boundaries,  the  concept  of  an 
array  of  dislocations  is  not  applicable,  and  there 
must  be  an  interface  across  which  foe  bonding 
departs  greatly  from  that  in  a  perfect  crystal.  For 
some  particular  orientations,  foe  bonding  may  be 
less  irregular  than  for  others,  and  the  grain  bound¬ 
ary  energy  will  depend  on  foe  relative  misorienta- 
tion  of  foe  grains  and  on  foe  location  of  the  boimd- 
ary  between  foem,  as  hasbeen  observed  by  Suzuki 
and  Kuroiwa  (1972). 

A  commonly  assumed  conditicm  for  foe  forma¬ 
tion  of  any  spedal  type  of  low  energy  boimdary  is 
foat  foe  lattice  points  for  foe  two  half-crystals 
should  match  up  with  one  anofoer  in  a  periodic 
way  along  foe  plane  of  foe  boundary.  The  Coinci¬ 
dence-Site  Lattice  (CSL)  model,  which  has  been 
described  in  relation  to  ice  by  Higashi  (1978)  and 
Hondoh  (1988),  has  this  property,  but  other  re¬ 
quirements  imposed  by  foe  model  are  not  general¬ 
ly  foou^t  to  be  significant  (Sutton  1984).  A  favor¬ 
able  configuration  in  ice  is  one  in  which  foe  grains 
are  rotated  relative  to  one  anofoer  by  34.1°  about 
foe  [1010]  direction,  and  such  a  boundary  in  which 
the  grains  are  joined  across  foeir  (1211)  planes  is 
illustrated  in  Figure22.HoiulohandHigashi  (1978) 
have  grown  bicrystals  containing  this  type  of 
boundary  and  have  shown  by  X-ray  topography 
foat  foe  boimdary  is  not  flat,  but  is  made  up  of 
facets  foat  make  small  angles  to  one  anofoer.  They 
propos^that  a  fecet  that  makes  a  small  angle  to  foe 
exact  ( 1 21 1 )  plane  will  contain  an  array  of  intrinsic 
grain  boimdary  dislcxations  of  foe  small  Burgers 
vector  foown  in  Figure  23;  between  these  disloca¬ 
tions  foe  lattices  will  match  exactly  as  in  Figure  22. 
Experiments  on  foe  diffusional  motion  of  boimd- 
aries  during  strain-fiee  annealing  have  shown 
variations  with  foe  type  of  boundary  and  some 
tendency  to  form  facets  wifoparticularorientations 
(Hondoh  and  Higashi  1979,  Nasello  et  al.  1992). 

Ofoer  examples  of  misorientations  foat  satisfy 
the  CSL  ccmditions  ate  47°  about  [1010]  and  21.8° 
about  [0001] .  The  latter  is  foought  to  lead  to  certain 
types  of  12-branched  snow  flakes  (Kobayashi  and 
^rukawa  1975).  These  special  cases  can  also  be 
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(00011,^^ 


Figun  22.  Boundary  on  ®  “o  ®  *o  ®  *o  “  "o 
a  (1211)  plane  between 

tzvo  grains  rotated  rela-  [OOQiY^^^ 

tmetooneanotherby34.1  ° 

(Aout  the  [lOTO]  axis  perpendicular  to  the  diagram.  The 
circles  are  latt^  points  not  molecules. 

thought  of  as  growtf  \  twins,  thou^  in  ice  they  are 
not  usually  formed  in  the  ways  commonly  associat¬ 
ed  with  twinning.  It  must  be  remember^  that  the 
coincidence  of  a  number  of  lattice  points  across  an 
interface  does  not  ensiue  foat  hydrogen  bonds  link 
up  in  a  manner  even  approximating  that  in  the 
perfect  crystal,  and  there  must  be  considerable  dis¬ 
order  of  foe  molecules  in  this  region.  The  thickness 
of  this  nonoystalline  or  'liquid-like'  region  is  not 
known.  In  this  connection  it  is  relevant  to  note  that, 
close  to  foe  melting  point,  liquid  water  is  present  in 
polycrystalline  ice  in  veins  that  lie  along  foe  lines  of 
intersection  of  grain  boimdaries  (Mader  1992,  Nye 
1992),  and  this  will  be  in  equilibrium  wifo  foe 
intenud  structure  of  foe  boundaries  themselves. 

Grain  boimdaries  can  moveinamunber  of  ways, 
according  to  principles  that  are  generally  appli^- 
ble  to  all  materials.  Under  a  shear  stress  a  pair  of 
suitably  oriented  grains  nuy  slide  over  one  another 
(HI  the  boundary  (Ignat  and  Frost 1987),  but  in  otlwr 
cases  shear  can  (xxrur  by  foe  migration  of  the  bound¬ 
ary  perpendict^  to  its  plane.  This  has  been  report- 
ed  for  a  34®  [1010]  boundary  by  Hondoh  (1988),  and 
explained  in  tennsof  the  glide  of  intrinsicboundary 
didocaticms  of  the  type  foown  in  Figure  23. 

The  moti(Hi  of  grain  boundaries,  or  of  dislrxsi- 
ticms  witiiin  them,  leads  to  a  characteristic  peak  in 
tile  low-frequency  internal  fricticm  (Perez  et  al. 
1979,  Tatibouet  et  al.  1987).  Bcnmdaries  also  move 
by  diffusive  prcxxsses  in  which  molecules  in  one 
grain  are  rearranged  in  foe  lattice  structure  of  foe 
otiier,  and  tills  is  what  happens  when  ice  recrystal¬ 
lizes  during  large-scale  plastic  flow. 

Being  places  where  tiiete  are  irregularities  in  foe 
lattke,  grain  boundaries  can  act  as  sources  and 
sinks  of  point  defects  and  dislocaticms.  We  have 
already  referred  to  the  formation  of  foulteddislcxa- 
tuxi  IcK^  by  the  ccmdensation  of  interstitials  on 
cooling.  Near  a  boundary  tiiete  is  a  zone  where 
these  kx^  are  not  formed  because  tile  excess  inter¬ 
stitials  are  lost  to  tile  bcHindary  (Liu  et  aL,  in  press). 
In  die  very  eariy  stages  of  deformaticHi,  stress  ccm- 


®  O**  ®  o  °  O  °  O  Figure  23.  Edge  dislo- 

*o  °  *oT^t  ^  •*  cation  in  the  boundary 

o  o  o  o  ^  produced  by  shear  in 

o  o  ®^  o  ®^  o  theplaneofthebound- 
ary  by  the  small  Burg¬ 
ers  vector  shown.  This  results  in  the  displacement  of  the 
boundary  downwards  in  the  right-hand  part  of  the 
diagram. 

centrations  form  at  grain  boundaries,  and  disloca¬ 
tion  IcKips  are  nucleated  from  them  (Hondoh  and 
Higashi  1983,  Liu  et  al.  1993). 

Electrical  properties  of  grain  boundaries 
in  doped  ice 

Electrical  measurements  on  monociystalline 
and  polycrystalline  samples  of  ice  lead  to  foe  con¬ 
clusion  that  these  types  of  ice  differ  significantly 
even  when  foe  concentrations  of  impurities  aver¬ 
aged  over  foe  volume  are  foe  same.  As  i(c  does  not 
exhibit  any  anisotropy  of  conductivity,  we  are  led 
to  conclude  that  foe  grain  boundaries  make  an 
appreciable  contribution  to  foe  conductivity. 

The  most  obvious  reason  why  this  may  be  so  is 
that  foe  impurities  segregate  to  the  boundaries, 
and  particularly  to  foe  triple  intersecti(Hi  lines 
where  tiiree  grain  boundaries  meet.  Mulvaney  et 
al.  (1988)  and  Wolff  et  al.  (1988)  investigated  foe 
impurity  distribution  in  polycrystalline  Antarctic 
ice  using  a  scanning  electron  microscope  wifo  a 
cold  stage  and  anX-raymicroanalysis  facility  wifo 
a  spatial  resoluticm  of  10  nm  and  a  detecticm  limit 
of  5  mM  (equivalent  to  490  ppm  for  H2SO4).  They 
found  that,  although  foe  volume  average  concen- 
traticHis  were  quite  small  (182  ppb  for  Na*^, 320ppb 
for  Q",  764  ppb  for  SO4  and  41  ppb  for  N(^ 
according  to  i(m-chromatography  and  atomic 
absoiptionspectroscopy),tiiec(HicentrationofSO^ 
at  triple  juncticms  was  25  M  witiiin  an  area  of  1  m^. 
This  concentraticHi  is  close  to  foe  eutectic  temper¬ 
ature  (4.9  M  freezing  at  -73®C),  so  tiiat  sudi  triple 
iuncti(His  remain  liquid  at  very  low  temperatures 
and  formanetworkof  filaments  of  extremely  hig^ 
electrical  conductiviQr.  Wolff  and  Paten  (1984) 
have  suggested  tiiat  foe  d.c.  ccmductivity  of  polar 
ice  could  be  caused  by  foe  presence  of  acidic  liquid 
layers  at  the  grain  boundaries.  They  showed  that 
it  is  plausible  that  these  impurities  will  concen¬ 
trate  at  the  triple  junctions,  and,  using  reliable  data 
for  H2SO4,  and  HCl  c(Hicentrati(His  at  foe 
Soufo  Pole,  they  derived  foe  correct  magnitude 
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and  temperature  dependence  for  the  conductivity 
ofsudiice. 

Ihe  di^rencea  between  dte  electrical  proper- 
tiesofmonoaystalline  and  polycrystalline  ke  seem 
to  persist  in  nominally  pure  samples.  This  sug¬ 
gests  that  die  grain  botmdaiies  may  have  an  en¬ 
hanced  conductivity  because  of  die  presence  of  a 
Tiquid-like' layer,  such  as  has  also  b^  proposed 
for  die  hree  surfece  (see  the  future  report  on  die 
Surface  of  Ice).  However,  at  present  this  is  just 
speculatiiHi,  and  experiments  (m  genuinely  pure 
boundaries  wiU  be  difficult  to  achieve. 
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